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Abstract 

We use the stability conditions of the Landau parameters for the symmetric nuclear matter and 
pure neutron matter to calculate the critical densities for the Skyrme type effective nucleon-nucleon 
interactions. We find that the critical density can be maximized by adjusting appropriately the 
values of the enhancement factor k associated with isovector giant dipole resonance, the quantity 
L which is directly related to the slope of the symmetry energy and the Landau parameter G' . 
However, restricting k, L and G' Q to vary within acceptable limits reduces the maximum value for 
the critical density p cr by ~ 25%. We also show that among the various quantities characterizing the 
symmetric nuclear matter, p cr depends strongly on the isoscalar effective mass m*/m and surface 
energy coefficient E s . For realistic values of m*/m and E s we get p cr = 2po to 3po (po = 0.16fm~ 3 ). 

PACS numbers: 21.30.Fe, 21.60.Jz,21.65.+f, 21.60. +c 
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There are many different parameterizations of the Skyrme type effective nucleon-nucleon 
interaction as recently reviewed in Ref. lj. Values of the Skyrme parameters are usually 
obtained by fitting the results of the Hartree-Fock calculations for the binding energies, 
charge radii and spin-orbit splittings of a few closed shell nuclei to the observed ones using 
a least square procedure. The Skyrme parameters are also constrained to yield appropriate 
values for some properties of the infinite symmetric nuclear matter at the saturation density 
Pnm- Over time, several improvements were made by appropriately treating the center of 
mass correction to the binding energy, modifying the spin-orbit part of the interaction and 
including in the least square fit the equation of state (EOS) for the pure neutron matter 
obtained from a realistic interaction DD . Nevertheless, not all the Skyrme parameters are 
well determined in this way. 

Recently, it has been shown in Ref. jj| that a more realistic parameterization of the 
Skyrme interaction can be obtained by subjecting it to stability requirements of the EOS 
defined by the inequality conditions for the Landau parameters for symmetric nuclear matter 
and pure neutron matter. In other words, few of the Skyrme parameters not well determined 
by the experimental data used in the least square procedure can be restricted by requiring 
that the inequality conditions are satisfied up to a maximum value for a nuclear matter 
density, also referred to as the critical density p cr . A very recent systematic study carried 
out in Ref. 3 using several Skyrme interactions indicates that the density dependence 
of the symmetry energy coefficient J plays a critical role in determining the properties of 
neutron star. Out of 87 different parameterizations for the Skyrme interaction considered 
in Ref. p only 27 of them, having a positive slope for the symmetry energy coefficient 
at nuclear matter densities p up to 3p (Po — 0.16fm~ 3 ), are found to be suitable for the 
neutron star model. Thus, it appears that the parameters of the Skyrme interactions can 
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be better constrained by combining the findings of Refs. 

In the present work we study the dependence of p cr on the nuclear matter saturation den- 
sity p nm , binding energy per nucleon B/A, isoscalar effective mass m*/m, incompressibility 
coefficient K nm , surface energy E s , and symmetry energy coefficient J, as also considered 
in Ref. 3- ft m ust be emphasized that we determine p cr in terms of the enhancement 
factor k, the coefficient L = 3pdJ/dp and the Landau parameter G' (at p nm ) instead of the 
combinations tjXj (i = 1,2, and 3) of the Skyrme parameters tj and Xi as used in Ref. "\. 
Unlike the combinations Uxi, the quantities k, L and G' Q , which can be expressed in terms 



2 



of the Skyrme parameters, are related to some physical processes. The enhancement factor 
k accounts for the deviations from the Thomas-Reiche-Kuhn (TRK) sum rule in the case 
of the isovector giant dipole resonance. The value of k at the p nm is expected to be ~ 0.5 
P, (|. The slope of the symmetry energy coefficient must be positive for p up to 3po, so 
that resulting Skyrme interaction can be suitable for the study of the properties of neutron 
stars p. The Landau parameter G' must be positive at p < p nm in order to appropriately 

n 

describe the position of the isovector Ml and Gamow- Teller states We find that by 
restricting the values of k, L and G' within the acceptable limits, the maximum critical 
densities are lowered by about 25% compared to the ones obtained without such restrictions 
^J. The constrains proposed in the present work not only maximizes the value of the critical 
density, but also ensures that the resulting Skyrme interaction can be used to study the bulk 
properties of finite nuclei as well as those of neutron stars. 

Recently a generalized form for the Skyrme energy density functional has been 

obtained using the Hohenberg-Kohn-Sham approach. However, in the present work we 
restrict ourselves to the energy density functional associated with the commonly used Skyrme 
interaction 0, 0] , 
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where, P° 2 is the spin exchange operator, a \ is the Pauli spin operator, k\ 2 = — z(Vi — 
V 2 )/2 and k \ 2 = — i(Vi — V 2 )/2. Here, the right and left arrows indicate that the mo- 
mentum operators act on the right and on the left, respectively. The Skyrme parameters ti, 
Xi and a for a fixed value of Wq can be expressed in terms of the quantities associated with 
the symmetric nuclear matter as follows Q> U, 0] . 
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In Eqs. (l2*|)- (jl(J|) . the various quantities characterizing the nuclear matter are the binding 
energy per nucleon B/A, isoscalar effective mass m*/m, nuclear matter incompressibility 
coefficient K nm , symmetry energy coefficient J, the coefficient L which is directly related 
to the slope of the symmetry energy coefficient (L = 3pdJ/dp), enhancement factor k and 
Landau parameter G' . All these quantities are taken at the saturation density p nm . It must 



be pointed out that we include consistently all terms in the energy density functional. In 
particular, the expression for the parameter G' used in the Eq. (@J is consistent with the 
energy density functional obtained using Eq. (fT|). For a more generalized Skyrme energy 
density functional, the expression for G' Q should be appropriately modified following Ref. 
7]. Once, T is known, T s can be calculated for a given value of the surface energy E s as 
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where, B(p)/A is the binding energy per nucleon given by, 
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The manner in which Eqs. (J2J) - (fTUj) can be used to evaluate the Skyrme parameters £», 
Xi and a is as follows. First, the parameters to and a can be obtained from Eqs. (0) and 
(JTUJ), respectively. Then, the parameter t 3 can be calculated using Eq. (jSJ). Next, T and 
T s can be calculated using Eqs. (|llj) and (|T4^) . respectively. Once, the combinations T and 
T s of the Skyrme parameters are known, one can calculate the remaining parameters in the 
following sequence, ti, t 2 , x 2 , Xi, x 3 and x . 
The stability criteria requires that [10], 
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where, X\ stands for the Landau parameters Fi, F(, Gi and G\ for a given multipolarity I. 
The Skyrme interaction only contains monopolar and dipolar contributions to the particle- 
hole interaction so that all Landau parameters are zero for I > 1. Thus, there are 12 different 
Landau parameters, i.e., F h F(, Gi and G\ (I = 0, 1) for the symmetric nuclear matter and 
F/ n) , G\ n) (I = 0,1) for the pure neutron matter. Each of these Landau parameters must 
satisfy the inequality condition given by Eq. ()17j) . Using the expressions for the Landau 
parameters in terms of the Skyrme parameters, given in Refs. |4| , one can obtain the values 
of the Landau parameters at any density for a given set of the Skyrme parameters. Thus, 
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the critical density which is nothing but the maximum density up to which all the inequality 
conditions are met can be easily determined. 

We have compiled the values of k, L and G' Q for several parameterization of the Skyrme 
interaction presented in Refs. [2,3,11-16]. We find that the values of k, L and G' Q vary 
over a wide ranges — 2, 40 — 160 MeV and —0.15 — 1.0, respectively. This is due to the 
fact that the experimental data used in the least-square procedure to fit the values of the 
Skyrme parameters can not constrain well the values of these quantities. These quantities 
can be constrained by requiring a reasonable value for the critical density. In what follows, 
we shall refer to the set of standard values for six quantities p nm = 0.16fm~ 3 , B/A = 16 
MeV, K nm = 230 MeV, m*/m = 0.7, E s = 18 MeV and J = 32MeV in short as the STD 
values for the nuclear matter input as used in Ref. jj]. 

We now study the dependence of p cr on the p nm , B/A, m*/m, K nm , E s and J. For a 
given set of values for these quantities, we calculate the maximum value of p cr by varying 
k, L and G' at the p nm within acceptable limits. We denote the maximum value of the 
critical density by p cr . It must be emphasized here that the present work differs from the 
one performed in Ref. by the fact that we constrain the values of k, L and G' , whereas, 
in Ref. jj], the value of p cr was calculated by varying the combinations tiXi (i = 1,2, and 
3) of Skyrme parameters with no restrictions. We find that if the space of the parameters 
tiXi is not restricted, it may lead to an unreasonable values of k and L. As an example, for 
the STD values of nuclear matter input, p cr becomes 3.5p for k = 1.0, L = 36 MeV and 
G'r, = 0.20 (at on). The value of G' n seems reasonable, but the value of k = 1.0 is a little too 

nn 

large |2J,|6J]. Further, we find that for p > p the value of L decreases with increasing p and 
it becomes negative for p > 1.6po? which makes the interaction not favorable for the neutron 
star model. 

In Fig. [T] we have displayed the results for p cr obtained by varying the various quantities 
associated with the nuclear matter around their standard values. To calculate p cr we allow 
for 0.25 < k < 0.5, < L < 100 MeV and < G' < 0.5 at the p nm . We further demand 
that L > at 3po- Fig. Q shows that p cr depends strongly on m*/m and E s . Whereas, p cr 
depends weakly on p nm , B/A and K nm and it is almost independent of J. These features for 
p cr are qualitatively similar to the ones presented in Ref. (4 1 . However, due to the restrictions 
imposed on the values of k, L and G' , the values of p cr becomes smaller than that obtained 
in Ref. j^J by up to 25%. For m*/m = 0.6 (0.7) and keeping all the other nuclear matter 
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quantities equal to their standard values, we get p cr = 4.5po (2.8po) compared to 6po (3.5po) 
obtained in Ref. 0]. 

It may be instructive to present the values of k, L and G' required to obtained p cr for a 
given set of values for the nuclear matter input. We find that k lies in the range of 0.45 — 0.5 
for variations in the nuclear matter input by up to ±15% relative to their standard values. 
These values of k clearly reflects the fact that restricting k to take values in the range of 
0.25 — 0.5 delimits the p cr to a lower value. In Figs. |2t and|2t> we have plotted the values 
of L and G' (at p nm ), respectively, which are needed to yield the p cr . We see that L varies 
from 20 to 60 MeV for different values of the nuclear matter input. For the STD values 
of nuclear matter input we find that L = 47 MeV at p = po- This value is quite large 
compared to the values of L — 35, 27 and 16 MeV associated with the the SkzO, Skzl and 
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Skz2 interactions |4j, respectively, which were obtained for the same standard values of the 
nuclear matter input, but by varying the combinations txx% and £2^2 with no restrictions 
and keeping £3X3 fixed to some arbitrary values. We see from Fig. EJd that, except for J, 
the value of G' Q (at p„ m ) depends strongly on the values of the various quantities associated 
with nuclear matter. The dependence of G' Q on the surface energy coefficient E s is the most 
pronounced one. We note that E s is mainly determined by the ground state properties of 
light nuclei. Thus, the center of mass correction to the binding energy and charge radii 
must be appropriately taken into account as they are very important for light nuclei and 
may affect the values obtained for E s . We see from Fig. |2)d that G' tends to vanish rapidly 
with increasing E s . 

As pointed out earlier (see from Fig. HJ), the dependence of p cr on p nm , J, and B/A is 
quite weak. Thus, it may be sufficient to calculate Per cLS cL function of m*/m, E s and K nm 
only. In Fig. |3] we display our results for the variation of E s versus m*/m, obtained for 
fixed values of p cr with the remaining nuclear matter quantities kept equal to their standard 
values. It can be seen from Fig. Elthat for a fixed value of p cr , E s decreases with the increase 
in m* lm. It is quite interesting to note that for E s — 18 ± 1 MeV (as most of the Skyrme 
interactions yield), p cr = 2po and 3po for m*/m = 0.72 — 0.85 and 0.63 — 0.73, respectively. 
For p cr = 4p one must have m* jm ~ 0.65 for not too low value of E s . The value of m* lm 



is also constrained by the centroid energy of the isoscalar giant quadrupole resonance [17 1 
which favors m* jm > 0.7. Thus, for reasonable values of E s and m*/m, one may obtain a 
Skyrme interaction with p cr = 2po to 3po- 
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In summary, we have used the stability conditions of the Landau parameters for the 
symmetric nuclear matter and pure neutron matter to calculate the critical densities for the 
Skyrme type effective nucleon-nucleon interactions. We find that the critical density can be 
maximized by appropriately adjusting the values of the enhancement factor k, coefficient L 
and the Landau parameter G' Q as these quantities are not well determined by the Skyrme 
parameters, conventionally obtained by fitting the experimental data for the ground state 
properties of finite nuclei. We exploit the fact that i) The value of k should be in the range 
of 0.25 — 0.5, needed to describe the TRK sum rule for the isovector giant dipole resonance 
, ii) L > for < p < 3po! a condition necessary for a Skyrme interaction to be suitable 
for studying the properties of neutron star Q]. and iii) G' > in order to reproduce the 
energies of the isovector Ml and Gamow- Teller states |6( . The maximum value of the critical 
density so obtained is lower by up to 25% compared to the ones obtained without any such 
restrictions |4{ . We show that the critical density obtained for realistic values of the surface 
energy coefficient (E s = 18 ± 1 MeV) and isoscalar effective mass (m*/m = 0.7 ± 0.1) lie in 
the range of 2po to 3po- Finally, we would like to remark that the pairing correlations were 
not included in the present work. It may be worthwhile to investigate the influence of the 
pairing correlations on the various quantities associated with the nuclear matter which may 
alter the value of the critical density. 
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FIG. 1: The dependence of critical density 

Per on the relative variation of p nm (dotted), B/A 
(dashed), m*/m (solid), K nm (open circles), E s (dashed-dot), and J (dashed- filled squares) around 
their standard values. 

FIG. 2: Variation of (a) the coefficient L and (b) the Landau parameter G' as a function of the 
various quantities associated with nuclear matter at p n m- The values of L and G' are determined 
by maximizing the critical density for a given set of values for the nuclear matter quantities. 

FIG. 3: Variations of the surface energy coefficient E s at p nm as a function of the effective mass 
m* /m for fixed values of the critical density p cr = 2po, 3po and 4po as labeled. All the other nuclear 
matter quantities are kept equals to their standard values. 
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